Genetic diversity of endosymbiotic bacteria Wolbachia infecting two mosquito species of the genus Eretmapodites occurring in sympatry in the Comoros archipelago.

Introduction: The influence of Wolbachia on mosquito reproduction and vector competence has led to renewed interest in studying the genetic diversity of these bacteria and the phenotypes they induced in mosquito vectors. In this study, we focused on two species of Eretmapodites, namely Eretmapodites quinquevittatus and Eretmapodites subsimplicipes, from three islands in the Comoros archipelago (in the Southwestern Indian Ocean). Methods: Using the COI gene, we examined the mitochondrial genetic diversity of 879 Eretmapodites individuals from 54 sites. Additionally, we investigated the presence and genetic diversity of Wolbachia using the wsp marker and the diversity of five housekeeping genes commonly used for genotyping through Multiple Locus Sequence Typing (MLST). Results and discussion: Overall, Er. quinquevittatus was the most abundant species in the three surveyed islands and both mosquito species occurred in sympatry in most of the investigated sites. We detected a higher mitochondrial genetic diversity in Er. quinquevittatus with 35 reported haplotypes (N = 615 specimens, Hd = 0.481 and π = 0.002) while 13 haplotypes were found in Er. subsimplicipes (N = 205 specimens, Hd = 0.338 and π = 0.001), this difference is likely due to the bias in sampling size between the two species. We report for the first time the presence of Wolbachia in these two Eretmapodites species. The prevalence of Wolbachia infection varied significantly between species, with a low prevalence recorded in Er. quinquevittatus (0.8%, N = 5/627) while infection was close to fixation in Er. subsimplicipes (87.7%, N = 221/252). Both male and female individuals of the two mosquito species appeared to be infected. The analysis of MLST genes revealed the presence of two Wolbachia strains corresponding to two new strain types (STs) within the supergroups A and B, which have been named wEretA and wEretB. These strains were found as mono-infections and are closely related, phylogenetically, to Wolbachia strains previously reported in Drosophila species. Finally, we demonstrate that maternal transmission of Wolbachia is imperfect in Er. subsimplicipes, which could explain the presence of a minority of uninfected individuals in the field.

Genetic characterization of dengue virus serotype 1 circulating in Reunion Island, 2019-2021, and the Seychelles, 2015-2016.

BACKGROUND: An unprecedent increase in the number of cases and deaths reported from dengue virus (DENV) infection has occurred in the southwestern Indian ocean in recent years. From 2017 to mid-2021 more than 70,000 confirmed dengue cases were reported in Reunion Island, and 1967 cases were recorded in the Seychelles from 2015 to 2016. Both these outbreaks displayed similar trends, with the initial circulation of DENV-2 which was replaced by DENV-1. Here, we aim to determine the origin of the DENV-1 epidemic strains and to explore their genetic characteristics along the uninterrupted circulation, particularly in Reunion. METHODS: Nucleic acids were extracted from blood samples collected from dengue positive patients; DENV-1 was identified by RT-qPCR. Positive samples were used to infect VERO cells. Genome sequences were obtained from either blood samples or infected-cell supernatants through a combination of both Illumina or MinION technologies. RESULTS: Phylogenetic analyses of partial or whole genome sequences revealed that all DENV-1 sequences from Reunion formed a monophyletic cluster that belonged to genotype I and were closely related to one isolate from Sri Lanka (OL752439.1, 2020). Sequences from the Seychelles belonged to the same major phylogenetic branch of genotype V, but fell into two paraphyletic clusters, with greatest similarity for one cluster to 2016-2017 isolate from Bangladesh, Singapore and China, and for the other cluster to ancestral isolates from Singapore, dating back to 2012. Compared to publicly available DENV-1 genotype I sequences, fifteen non-synonymous mutations were identified in the Reunion strains, including one in the capsid and the others in nonstructural proteins (NS) (three in NS1, two in NS2B, one in NS3, one in NS4B, and seven in NS5). CONCLUSION: In contrast to what was seen in previous outbreaks, recent DENV-1 outbreaks in Reunion and the Seychelles were caused by distinct genotypes, all likely originating from Asia where dengue is (hyper)endemic in many countries. Epidemic DENV-1 strains from Reunion harbored specific non-synonymous mutations whose biological significance needs to be further investigated.

Evidence of Eretmapodites subsimplicipes and Aedes albopictus as competent vectors for Rift Valley fever virus transmission in Mayotte.

Rift Valley fever (RVF) recently re-emerged in Mayotte. We described, for the first time, that the mosquito species Eretmapodites subsimplicipes, a highly abundant species in Mayotte, is a competent vector for the transmission of RVF virus using three distinct populations native to Mayotte. We also showed that Aedes albopictus specimens are able to transmit RVF virus (RVFV) as previously observed in mosquito populations of other countries emphasizing the need of the increase vigilance for this highly invasive species of global distribution. Altogether, these results underline the epidemiological importance of both species for RVFV transmission in Mayotte and contribute to better understand the RVF epidemiological cycle and help to implement efficient prevention measures.

Fitness costs associated with a GABA receptor mutation conferring dieldrin resistance in Aedes albopictus.

Understanding the dynamics of insecticide resistance genes in mosquito populations is pivotal for a sustainable use of insecticides. Dieldrin resistance in Aedes albopictus is conferred by the alanine to serine substitution (A302S or RdlR allele) in the γ-aminobutyric acid (GABA) receptor encoded by the Rdl gene. On Reunion Island, dieldrin resistance was initially reported in natural Ae. albopictus populations sampled in 2008 despite the ban of dieldrin since 1994. To monitor insecticide resistance in Ae. albopictus on the island and to identify its drivers, we measured (i) the frequency of resistance alleles in 19 distinct natural populations collected between 2016 and 2017, (ii) fitness costs associated with dieldrin resistance in laboratory-controlled experiments, and (iii) the resistance conferred by RdlR to fipronil, an insecticide widely used on the island and reported to cross-react with RdlR. The results show a persistence of RdlR in Ae. albopictus natural populations at low frequencies. Among the measured life history traits, mortality in pre-imaginal stages, adults' survival as well as the proportion of egg-laying females were significantly affected in resistant mosquitoes. Finally, bioassays revealed resistance of RdlR mosquitoes to fipronil, suggesting that the use of fipronil in natura could select for the RdlR allele. This study shows that dieldrin resistance is persistent in natural mosquito populations likely as a result of combined effects between fitness costs associated with RdlR and selection exerted by cross-reacting environmental insecticides such as fipronil.

Zika vector competence data reveals risks of outbreaks: the contribution of the European ZIKAlliance project.

First identified in 1947, Zika virus took roughly 70 years to cause a pandemic unusually associated with virus-induced brain damage in newborns. Zika virus is transmitted by mosquitoes, mainly Aedes aegypti, and secondarily, Aedes albopictus, both colonizing a large strip encompassing tropical and temperate regions. As part of the international project ZIKAlliance initiated in 2016, 50 mosquito populations from six species collected in 12 countries were experimentally infected with different Zika viruses. Here, we show that Ae. aegypti is mainly responsible for Zika virus transmission having the highest susceptibility to viral infections. Other species play a secondary role in transmission while Culex mosquitoes are largely non-susceptible. Zika strain is expected to significantly modulate transmission efficiency with African strains being more likely to cause an outbreak. As the distribution of Ae. aegypti will doubtless expand with climate change and without new marketed vaccines, all the ingredients are in place to relive a new pandemic of Zika.

Overview of dengue outbreaks in the southwestern Indian Ocean and analysis of factors involved in the shift toward endemicity in Reunion Island: A systematic review.

BACKGROUND: Dengue is the world's most prevalent mosquito-borne viral disease. It is endemic in many tropical and subtropical countries and represents a significant global health burden. The first reports of dengue virus (DENV) circulation in the South West Indian Ocean (SWIO) islands date back to the early 1940s; however, an increase in DENV circulation has been reported in the SWIO in recent years. The aim of this review is to trace the history of DENV in the SWIO islands using available records from the Comoros, Madagascar, Mauritius, Mayotte, Seychelles, and Reunion. We focus in particular on the most extensive data from Reunion Island, highlighting factors that may explain the observed increasing incidence, and the potential shift from one-off outbreaks to endemic dengue transmission. METHODS: Following the PRISMA guidelines, the literature review focused queried different databases using the keywords "dengue" or "Aedes albopictus" combined with each of the following SWIO islands the Comoros, Madagascar, Mauritius, Mayotte, Seychelles, and Reunion. We also compiled case report data for dengue in Mayotte and Reunion in collaboration with the regional public health agencies in these French territories. References and data were discarded when original sources were not identified. We examined reports of climatic, anthropogenic, and mosquito-related factors that may influence the maintenance of dengue transmission independently of case importation linked to travel. FINDINGS AND CONCLUSIONS: The first report of dengue circulation in the SWIO was documented in 1943 in the Comoros. Then not until an outbreak in 1976 to 1977 that affected approximately 80% of the population of the Seychelles. DENV was also reported in 1977 to 1978 in Reunion with an estimate of nearly 30% of the population infected. In the following 40-year period, DENV circulation was qualified as interepidemic with sporadic cases. However, in recent years, the region has experienced uninterrupted DENV transmission at elevated incidence. Since 2017, Reunion witnessed the cocirculation of 3 serotypes (DENV-1, DENV-2 and DENV-3) and an increased number of cases with severe forms and deaths. Reinforced molecular and serological identification of DENV serotypes and genotypes circulating in the SWIO as well as vector control strategies is necessary to protect exposed human populations and limit the spread of dengue.

Spatio-Temporal Dynamics of a Dieldrin Resistance Gene in Aedes albopictus and Culex quinquefasciatus Populations From Reunion Island.

The control of mosquito populations using insecticides is increasingly threatened by the spread of resistance mechanisms. Dieldrin resistance, conferred by point mutations in the Rdl gene encoding the γ-aminobutyric acid receptor, has been reported at high prevalence in mosquito populations in response to selective pressures. In this study, we monitored spatio-temporal dynamics of the resistance-conferring RdlR allele in Aedes (Stegomyia) albopictus (Skuse, 1895) and Culex (Culex) quinquefasciatus (Say, 1823) populations from Reunion Island. Specimens of both mosquito species were sampled over a 12-month period in three cities and in sites located at lower (<61 m) and higher (between 503 and 564 m) altitudes. Mosquitoes were genotyped using a molecular test detecting the alanine to serine substitution (A302S) in the Rdl gene. Overall, the RdlR frequencies were higher in Cx. quinquefasciatus than Ae. albopictus. For both mosquito species, the RdlR frequencies were significantly influenced by location and altitude with higher RdlR frequencies in the most urbanized areas and at lower altitudes. This study highlights environmental factors that influence the dynamics of insecticide resistance genes, which is critical for the management of insecticide resistance and the implementation of alternative and efficient vector control strategies.

Wolbachia modulates prevalence and viral load of Culex pipiens densoviruses in natural populations.

The inadequacy of standard mosquito control strategies calls for ecologically safe novel approaches, for example the use of biological agents such as the endosymbiotic α-proteobacteria Wolbachia or insect-specific viruses (ISVs). Understanding the ecological interactions between these "biocontrol endosymbionts" is thus a fundamental step. Wolbachia are transmitted vertically from mother to offspring and modify their hosts' phenotypes, including reproduction (e.g., cytoplasmic incompatibility) and survival (e.g., viral interference). In nature, Culex pipiens (sensu lato) mosquitoes are always found infected with genetically diverse Wolbachia called wPip that belong to five phylogenetic groups. In recent years, ISVs have also been discovered in these mosquito species, although their interactions with Wolbachia in nature are unknown. Here, we studied the interactions between a widely prevalent ISV, the Culex pipiens densovirus (CpDV, Densovirinae), and Wolbachia in northern Tunisian C. pipiens populations. We showed an influence of different Wolbachia groups on CpDV prevalence and a general positive correlation between Wolbachia and CpDV loads. By investigating the putative relationship between CpDV diversification and wPip groups in the different sites, we detected a signal linked to wPip groups in CpDV phylogeny in sites where all larvae were infected by the same wPip group. However, no such signal was detected where the wPip groups coexisted, suggesting CpDV horizontal transfer between hosts. Overall, our results provide good evidence for an ecological influence of Wolbachia on an ISV, CpDV, in natural populations and highlight the importance of integrating Wolbachia in our understanding of ISV ecology in nature.

Contrasted transmission efficiency of Zika virus strains by mosquito species Aedes aegypti, Aedes albopictus and Culex quinquefasciatus from Reunion Island.

BACKGROUND: Zika virus (ZIKV) is a mosquito-borne flavivirus that recently emerged in the South Pacific islands and Americas where unprecedented outbreaks were reported. Although Aedes aegypti is considered to be the main vector for ZIKV, other mosquito species have been shown to be potential vectors and differences in vector competence with respect to mosquito strain and ZIKV strain have been demonstrated. In this study we compared the vector competence of three mosquito species Aedes albopictus, Ae. aegypti and Culex quinquefasciatus from Reunion Island for three ZIKV strains. METHODS: Five mosquito strains (2 strains of Ae. albopictus, 1 of Ae. aegypti and 2 of Cx. quinquefasciatus) were exposed to three ZIKV strains: one African strain (Dak84) and two Asian strains (PaRi_2015 and MAS66). The vector competence parameters (infection rate, dissemination efficiency and transmission efficiency) and viral loads were examined at 14 and 21 days post-infection. RESULTS: The two Cx. quinquefasciatus strains did not become infected and were therefore unable to either disseminate or transmit any of the three ZIKV strains. Aedes albopictus and Ae. aegypti strains were poorly competent for the two Asian ZIKV strains, while both mosquito species displayed higher infection rates, dissemination and transmission efficiencies for the African ZIKV Dak84 strain. However, this African ZIKV strain was better transmitted by Ae. aegypti as compared to Ae. albopictus. CONCLUSIONS: Our results show that both Ae. albopictus and Ae. aegypti, from Reunion Island, are more likely to be competent for ZIKV in contrast to Cx. quinquefasciatus which appeared refractory to all tested ZIKV strains. This improves our understanding of the role of mosquito species in the risk of the ZIKV emergence on Reunion Island.

High genetic diversity but no geographical structure of Aedes albopictus populations in Réunion Island.

BACKGROUND: In recent years, the Asian tiger mosquito Aedes albopictus has emerged as a species of major medical concern following its global expansion and involvement in many arbovirus outbreaks. On Réunion Island, Ae. albopictus was responsible for a large chikungunya outbreak in 2005-2006 and more recently an epidemic of dengue which began at the end of 2017 and is still ongoing at the time of writing. This dengue epidemic has seen a high number of human cases in south and west coastal regions, while few cases have been reported in the north and east of the island. To better understand the role of mosquito populations in such spatial patterns of dengue virus transmission in Réunion Island, we examined the genetic diversity and population structure of Ae. albopictus sampled across the island. RESULTS: Between November 2016 and March 2017, a total of 564 mosquitoes were collected from 19 locations in three main climatic regions (West, East and Center) of Réunion Island and were genotyped using 16 microsatellite loci. A high genetic diversity was observed with 2-15 alleles per locus and the average number of alleles per population varying between 4.70-5.90. Almost all FIS values were significantly positive and correlated to individual relatedness within populations using a hierarchical clustering approach based on principal components analyses (HCPC). However, the largest part of genetic variance was among individuals within populations (97%) while only 3% of genetic variance was observed among populations within regions. Therefore, no distinguishable population structure or isolation by distance was evidenced, suggesting high rates of gene flow at the island scale. CONCLUSIONS: Our results show high genetic diversity but no genetic structure of Ae. albopictus populations in Réunion Island thus reflecting frequent movements of mosquitoes between populations probably due to human activity. These data should help in the understanding of Ae. albopictus vector capacity and the design of effective mosquito control strategies.

Insecticide resistance genes affect Culex quinquefasciatus vector competence for West Nile virus.

Insecticide resistance has been reported to impact the interactions between mosquitoes and the pathogens they transmit. However, the effect on vector competence for arboviruses still remained to be investigated. We examined the influence of two insecticide resistance mechanisms on vector competence of the mosquito Culex quinquefasciatus for two arboviruses, Rift Valley Fever virus (RVFV) and West Nile virus (WNV). Three Cx. quinquefasciatus lines sharing a common genetic background were used: two insecticide-resistant lines, one homozygous for amplification of the Ester2 locus (SA2), the other homozygous for the acetylcholinesterase ace-1 G119S mutation (SR) and the insecticide-susceptible reference line Slab. Statistical analyses revealed no significant effect of insecticide-resistant mechanisms on vector competence for RVFV. However, both insecticide resistance mechanisms significantly influenced the outcome of WNV infections by increasing the dissemination of WNV in the mosquito body, therefore leading to an increase in transmission efficiency by resistant mosquitoes. These results showed that insecticide resistance mechanisms enhanced vector competence for WNV and may have a significant impact on transmission dynamics of arboviruses. Our findings highlight the importance of understanding the impacts of insecticide resistance on the vectorial capacity parameters to assess the overall consequence on transmission.

Construction of a genetic sexing strain for Aedes albopictus: a promising tool for the development of sterilizing insect control strategies targeting the tiger mosquito.

BACKGROUND: Aedes albopictus is an invasive mosquito species of global medical concern as its distribution has recently expanded to Africa, the Americas and Europe. In the absence of prophylaxis protecting human populations from emerging arboviruses transmitted by this mosquito species, the most straightforward control measures rely on the suppression or manipulation of vector natural populations. A number of environmental-friendly methods using mass releases of sterilizing males are currently under development. However, these strategies are still lacking an efficient sexing method required for production of males at an industrial scale. RESULTS: We present the first Genetic Sexing Strain (GSS) in Ae. albopictus, hereafter referred as Tikok, obtained by sex linkage of the rdl gene conferring dieldrin resistance. Hatch rate, larval survival and sex ratio were followed during twelve generations. The use of dieldrin at the third larval stage allowed selecting 98 % of males on average. CONCLUSION: A good production rate of Tikok males makes this GSS suitable for any control method based on mass production of Ae. albopictus males. Despite limitations resulting from reduced egg hatch as well as the nature of the used insecticide, the construction of this GSS paves the way for industrial sex separation of Ae. albopictus.

Author Correction: Culex pipiens crossing type diversity is governed by an amplified and polymorphic operon of Wolbachia.

In the originally published HTML and PDF versions of this Article, gel images in Figures 7c and 8c were not prepared as per the Nature journal policy. These figure panels have now been corrected in both the PDF and HTML versions of the Article.In Fig. 7c, the lane labelled 'Ha' was inappropriately duplicated to represent the lane labelled 'Ich13'. The corrected version of Fig. 7c includes PCR-RFLP on DNA from the Ichkeul 13 line, which had been run on a separate gel. The original unprocessed gel images are provided in Supplementary Figure 1 associated with this correction, with the relevant corresponding bands denoted. A repeat experiment of the PCR-RFLP test is also presented as Supplementary Figure 2.In Fig. 8c, the image was assembled from two separate gels without clear demarcation. The corrected Fig. 8c clearly separates lanes from the two gels, and the original unprocessed gel images are provided in the Supplementary Information associated with this correction.These corrections do not alter the original meaning of the experiments, their results, their interpretation, or the conclusions of the paper. We apologize for any confusion this may have caused to the readers of Nature Communications.

Discovery of Culex pipiens associated tunisia virus: a new ssRNA(+) virus representing a new insect associated virus family.

In the global context of arboviral emergence, deep sequencing unlocks the discovery of new mosquito-borne viruses. Mosquitoes of the species Culex pipiens, C. torrentium, and C. hortensis were sampled from 22 locations worldwide for transcriptomic analyses. A virus discovery pipeline was used to analyze the dataset of 0.7 billion reads comprising 22 individual transcriptomes. Two closely related 6.8 kb viral genomes were identified in C. pipiens and named as Culex pipiens associated tunisia virus (CpATV) strains Ayed and Jedaida. The CpATV genome contained four ORFs. ORF1 possessed helicase and RNA-dependent RNA polymerase (RdRp) domains related to new viral sequences recently found mainly in dipterans. ORF2 and 4 contained a capsid protein domain showing strong homology with Virgaviridae plant viruses. ORF3 displayed similarities with eukaryotic Rhoptry domain and a merozoite surface protein (MSP7) domain only found in mosquito-transmitted Plasmodium, suggesting possible interactions between CpATV and vertebrate cells. Estimation of a strong purifying selection exerted on each ORFs and the presence of a polymorphism maintained in the coding region of ORF3 suggested that both CpATV sequences are genuine functional viruses. CpATV is part of an entirely new and highly diversified group of viruses recently found in insects, and that bears the genomic hallmarks of a new viral family.

Culex pipiens crossing type diversity is governed by an amplified and polymorphic operon of Wolbachia.

Culex pipiens mosquitoes are infected with Wolbachia (wPip) that cause an important diversity of cytoplasmic incompatibilities (CIs). Functional transgenic studies have implicated the cidA-cidB operon from wPip and its homolog in wMel in CI between infected Drosophila males and uninfected females. However, the genetic basis of the CI diversity induced by different Wolbachia strains was unknown. We show here that the remarkable diversity of CI in the C. pipiens complex is due to the presence, in all tested wPip genomes, of several copies of the cidA-cidB operon, which undergoes diversification through recombination events. In 183 isofemale lines of C. pipiens collected worldwide, specific variations of the cidA-cidB gene repertoires are found to match crossing types. The diversification of cidA-cidB is consistent with the hypothesis of a toxin-antitoxin system in which the gene cidB co-diversifies with the gene cidA, particularly in putative domains of reciprocal interactions.

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment, sublethality and complementation.

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace-1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace-1 alleles (single-copy, susceptible S and insecticide-resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing-overs between local single-copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long-term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.

Comparison of Irradiation and Wolbachia Based Approaches for Sterile-Male Strategies Targeting Aedes albopictus.

The global expansion of Aedes albopictus together with the absence of vaccines for most of the arboviruses transmitted by this mosquito has stimulated the development of sterile-male strategies aiming at controlling disease transmission through the suppression of natural vector populations. In this context, two environmentally friendly control strategies, namely the Sterile Insect Technique (SIT) and the Wolbachia-based Incompatible Insect Technique (IIT) are currently being developed in several laboratories worldwide. So far however, there is a lack of comparative assessment of these strategies under the same controlled conditions. Here, we compared the mating capacities, i.e. insemination capacity, sterilization capacity and mating competitiveness of irradiated (35 Gy) and incompatible Ae. albopictus males at different ages and ratios under laboratory controlled conditions. Our data show that there was no significant difference in insemination capacity of irradiated and incompatible males, both male types showing lower capacities than untreated males at 1 day but recovering full capacity within 5 days following emergence. Regarding mating competitiveness trials, a global observed trend is that incompatible males tend to induce a lower hatching rate than irradiated males in cage controlled confrontations. More specifically, incompatible males were found more competitive than irradiated males in 5:1 ratio regardless of age, while irradiated males were only found more competitive than incompatible males in the 1:1 ratio at 10 days old. Overall, under the tested conditions, IIT seemed to be slightly more effective than SIT. However, considering that a single strategy will likely not be adapted to all environments, our data stimulates the need for comparative assessments of distinct strategies in up-scaled conditions in order to identify the most suitable and safe sterilizing technology to be implemented in a specific environmental setting and to identify the parameters requiring fine tuning in order to reach optimal release conditions.

Wolbachia-based population control strategy targeting Culex quinquefasciatus mosquitoes proves efficient under semi-field conditions.

In mosquitoes, the maternally inherited bacterial Wolbachia induce a form of embryonic lethality called cytoplasmic incompatibility (CI). This property can be used to reduce the density of mosquito field populations through inundative releases of incompatible males in order to sterilize females (Incompatible Insect Technique, or IIT, strategy). We have previously constructed the LR[wPip(Is)] line representing a good candidate for controlling field populations of the Culex quinquefasciatus mosquito in the islands of the south-western Indian Ocean. The main purpose of the present study was to fill the gap between laboratory experiments and field implementation, i.e. assessing mating competitiveness of these incompatible males under semi-field conditions. In a first set of experiments, we analyzed crossing relationships between LR[wPip(Is)] males and La Réunion field females collected as larvae in 19 distinct localities throughout the island. This investigation revealed total embryonic mortality, confirming the strong sterilizing capacity of LR[wPip(Is)] males. Subsequently, mating competitiveness of LR[wPip(Is)] males was assessed under semi-field conditions in the presence of field males and females from La Réunion. Confrontations were carried out in April and December using different ratios of LR[wPip(Is)] to field males. The results indicated that the LR[wPip(Is)] males successfully compete with field males in mating with field females, displaying even higher competitiveness than field males in April. Our results support the implementation of small-scale field tests in order to assess the feasibility of IIT against Cx. quinquefasciatus in the islands of southwestern Indian Ocean where this mosquito species is a proven competent vector for human pathogens.

Stable coexistence of incompatible Wolbachia along a narrow contact zone in mosquito field populations.

In arthropods, the intracellular bacteria Wolbachia often induce cytoplasmic incompatibility (CI) between sperm and egg, which causes conditional embryonic death and promotes the spatial spread of Wolbachia infections into host populations. The ability of Wolbachia to spread in natural populations through CI has attracted attention for using these bacteria in vector-borne disease control. The dynamics of incompatible Wolbachia infections have been deeply investigated theoretically, whereas in natural populations, there are only few examples described, especially among incompatible infected hosts. Here, we have surveyed the distribution of two molecular Wolbachia strains (wPip11 and wPip31) infecting the mosquito Culex pipiens in Tunisia. We delineated a clear spatial structure of both infections, with a sharp contact zone separating their distribution areas. Crossing experiments with isofemale lines from different localities showed three crossing types: wPip11-infected males always sterilize wPip31-infected females; however, while most wPip31-infected males were compatible with wPip11-infected females, a few completely sterilize them. The wPip11 strain was thus expected to spread, but temporal dynamics over 7 years of monitoring shows the stability of the contact zone. We examined which factors may contribute to the observed stability, both theoretically and empirically. Population cage experiments, field samples and modelling did not support significant impacts of local adaptation or assortative mating on the stability of wPip infection structure. By contrast, low dispersal probability and metapopulation dynamics in the host Cx. pipiens probably play major roles. This study highlights the need of understanding CI dynamics in natural populations to design effective and sustainable Wolbachia-based control strategies.

Wolbachia divergence and the evolution of cytoplasmic incompatibility in Culex pipiens.

Many insect species harbor Wolbachia bacteria that induce cytoplasmic incompatibility (CI), i.e. embryonic lethality in crosses between infected males and uninfected females, or between males and females carrying incompatible Wolbachia strains. The molecular mechanism of CI remains unknown, but the available data are best interpreted under a modification-rescue model, where a mod function disables the reproductive success of infected males' sperm, unless the eggs are infected and express a compatible resc function. Here we examine the evolution of CI in the mosquito Culex pipiens, harbouring a large number of closely related Wolbachia strains structured in five distinct phylogenetic groups. Specifically, we used a worldwide sample of mosquito lines to assess the hypothesis that genetic divergence should correlate with the divergence of CI properties on a low evolutionary scale. We observed a significant association of Wolbachia genetic divergence with CI patterns. Most Wolbachia strains from the same group were compatible whereas those from different groups were often incompatible. Consistently, we found a strong association between Wolbachia groups and their mod-resc properties. Finally, lines from the same geographical area were rarely incompatible, confirming the conjecture that the spatial distribution of Wolbachia compatibility types should be constrained by selection. This study indicates a clear correlation between Wolbachia genotypes and CI properties, paving the way toward the identification of the molecular basis of CI through comparative genomics.